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Abstract 



Riemann conjectured that all the zeros of the Riemann S-function are real, which is now 
known as the Riemann Hypothesis (RH) . In this article we introduce the study of the zeros 
of the truncated sums Sjv(z) in Riemann's uniformly convergent infinite series expansion of 
S(z) involving incomplete gamma functions. We conjecture that when the zeros of Sjv(z) 
in the first quadrant of the complex plane are listed by increasing real part, their imaginary 
parts are monotone nondecreasing. We show how this conjecture implies the RH, and discuss 
some computational evidence for this and other related conjectures. 

1 Introduction 

Following Riemann (as described in a copy of an English translation of his memoir contained in 
the appendix of [EdwOl]), let 



3(z) is an even, entire function, and the famous Riemann Hypothesis (RH) says that all the 
zeros of 3 are real. Let Q denote the first quadrant of the complex plane $t(z) > 0, > 0. 
Since 2 (2) is even and real on the real line, we can restate the RH as saying all zeros of 2(z) in 
Q are real. Since £(s) is nonzero in 3?(s) > 1, it follows that all zeros of 3 in Q satisfy Q(z) < .5. 
In 1914 Hardy (as reprinted in [BCRW08J) showed that 3(z) has infinitely many real zeros and 
in 1942 Selberg [Sel42] showed that a positive proportion of the zeros of 3(z) are real. More 
recent work of Conrey [Con89] has at least 2/5 of the zeros on the real line. 
Riemann derived the following expression for S(z); 






(2) 



where </>(*) = Y^Li 4>n(t) with 



(f> n (t) = exp(-n 2 7rexp(2t))(87r 2 n 4 exp(4.5t) - 12vrn 2 exp(2.5i)). 



(3) 
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The function <p(t) is known to be an even function of t. Polya I '6120 investigated ways of 
approximating (j){t) by simpler functions. He showed that if in ([2]) we replace </>(i) by 



0i(£) = exp(-7rcosh(2t))(8vr 2 cosh(4.5i) - 12vr cosh(2.5t)) 



(4) 



(obtained by replacing most of the exponentials in the definition of (fti(t) by hyperbolic cosines), 
then the resulting integral has only real zeros. Polya also showed that if we replace <p(t) by any 
function which is not an even function of t, then the resulting integral has only finitely many real 
zeros. Hejhal [Hej90] investigated what happens if we replace (fr(t) by a "Polya approximate", 
i.e. a finite sum of the form 

N 



By building on earlier work of Bombieri and Hejhal [BH87] . which showed that, under the Gen- 
eralized Riemann Hypothesis (GRH) and some other assumptions, certain linear combinations 
of Dirichlet L-series asymptotically have 100% of their zeros on the critical line, Hejhal was able 
to show that if we replace (f>(t) in ([2]) by a Polya approximate the resulting function asymptoti- 
cally also has 100% of its zeros on the real line (but infinitely many zeros off the line). By 100% 
asymptotically we mean that the proportion of zeros in Q satisfying $t(z) < m that are on the 
real line approaches 1 as m —* oo. 

The starting point for this investigation is the idea that perhaps it is not necessary for 
worthwhile approximates to have all their zeros on the real line. If a given family of approximates 
approach S(z) uniformly, and if for each element in the family one could prove that within a 
certain sub-region of Q all the zeros are real, with the size of the sub-region expanding to 
eventually include all of Q as our approximates approach H, then this would also imply RH. 
Thus it may be worth studying replacements for (j){t) in ([2]) which are not even. With this in 
mind, a natural question to ask is what happens if we replace <p(t) by Xm=i ^nif)- 

2 Preliminary calculations 

Let G(z; a, b) denote the integral 



where z € C, a, b £ M with a > 0. Making the change of variable t = aexp(2u), so dt = 
aexp(2u)2du, and du = dt/2t, we get 




n 2 7rcosh(2f))(87r 2 n 4 cosh(4.5t) - 127m 2 cosh(2.5t)). 



(5) 



71=1 






dt 



2t 



(7) 
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a' 



(10) 



(8) 



(9) 
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where 

T(z,a) = J eM~t)t z j (11) 

is the (upper) incomplete gamma function. For lack of a better name, we will refer to G(z; a, b) 
as a "hyperbolic gamma function". 
From ([2]) and © we have 

oo 

E(z) = ^2 2ir 2 n 4 G{z/2; n 2 vr, 9/4) - tom 2 G(z/2; n 2 vr, 5/4), (12) 

71=1 

the interchange in integration and summation being justified by the uniform convergence. For 
a £ R, a > 0, the function r(z, a) is entire (as a function of 2), and hence so is G(z; a, b). There 
is a routine in Maple to compute T(z, a), and (in the RootFinding package) a routine to compute 
the zeros, using the argument principle and Newton's method, of a given analytic function in any 
rectangle of the complex plane. Using this, the author made several computations to compute 
the zeros of the E-approximates 

JV 

5jv(s) (13) 

71=1 

where 

:= 2vr 2 n 4 G(z/2; n 2 vr, 9/4) - 3vrn 2 G(z/2; n 2 vr, 5/4), (14) 

for various small values of N. Lists of zeros for some of these are contained in the Appendix. 
In these computer runs, the parameter "Digits" in Maple (which tells the computer to use this 
many significant digits in all calculations) was typically set to 20A — 10 or so, whatever number 
of digits was needed to compute the function in question over the specified rectangle accurately 
to 20 or so significant digits. After runs were made first with Digits equal to 20N — 10, they 
were sometimes run again with Digits equal to 20 A, and the resulting zeros typically agreed to 
at least 16 decimal digits or so, which the author has taken to mean the computer generated 
zeros (for Digits equal to 20A — 10) agree with the actual ones to at least 10 decimal digits, 
although no attempt has been made to establish rigorous error bounds. 

We say that a given function F(z) has monotonic zeros in a region D of the complex plane 
if, when we list the zeros of F in D by increasing real part, the imaginary parts of the zeros are 
monotone nondecreasing. Formally, if {a±,a2, • • •} are the zeros of F in D numbered so that 
M(oti) < 3f?(aj + i for i > 1, then 9(a«) < 9(aj+i). (We assume F has at most one zero on the 
intersection of any vertical line with D.) The data in the Appendix and other computer runs 
support the following hypothesis. 

Conjecture 1 For JVgN, Hjv(z) has monotonic zeros in Q. 
Proposition 1 Conjecture U\ implies the Riemann Hypothesis. 

Proof. This follows from the argument principle, combined with the simple fact that a function 
with infinitely many positive real zeros, and with monotonic zeros in Q, has only real zeros in 
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Q. Assume the RH is false, and let r be the zero of E(z) in Q with minimal real part, among 
those zeros with positive imaginary part, and let r = cr + it. By the argument principle, we have 

' 1 E ' {Z) dz = 1, (15) 



2iri J Ce E(z) 

where the integral is taken counterclockwise around a circle C e centered at t, of small radius e, 
so no other zeros of E{z) are enclosed in C. Next choose iV sufficiently large so that En(z) has 
a real zero 7 with 7 > 2a. We can do this since Hat(z) converges uniformly to E{z) on compacta 
in Q, both E(z) and En(z) are real on the real line, and since E(z) has infinitely many positive 
real zeros. By assumption En(z) has monotonic zeros in Q, hence has no non-real zeros in Q 
with real part less than 2a. This implies 

/ §44^ = o> (is) 

and so 

1 r E'(z) E' N (z 



2iri J Cf E(z) E N (z 
1 X E'(z)E N (z)-E' N (z)E(z) 



dz (17) 

r _ _ dz. (18) 

2itiJ C€ E(z)E N {z) v ; 

On the closed and bounded set C e , \E(z) \ is nonzero and hence must assume an absolute minimum 
5 > 0. Due to the uniform convergence, as N — > 00, the minimum of |Hjy(z)| on C e must 
eventually be greater than 5/2. Hence for large iV the modulus of the denominator of the 
integrand in fp~8|) is bounded away from zero, but (since E' N (z) approaches E'(z) uniformly) the 
numerator approaches zero, and so the integral will also approach zero, a contradiction. □ 

Remark 1 A weaker form of Conjecture^ which still implies RH, is that there are no non-real 
zeros ofE]\i(z) in Q whose real part is less than the largest real zero ofE]^(z). Since En{z) is 
real for real z, the real zeros can be found by looking at sign changes along the real line. Then the 
argument principle can be used via a numerical integration to obtain the total number of zeros 
o/Hjy(z) with real part not greater than the largest real zero, and matched against the number 
of real zeros. Computations along these lines indicate this weaker form of Conjecture U\ is true 
at least for N < 10. Below we list the largest real zero and number of real zeros of Ejy(z) for 
N < 10. 



K 


largest real zero 


number of real zeros 


1 


14.0454395788 


1 


2 


39.5324810798 


7 


3 


65.0320737720 


15 


4 


103.3679880094 


32 


5 


149.0026994921 


53 


6 


197.9575955732 


79 


7 


258.5304836632 


113 


8 


327.3794646017 


155 


9 


406.8174206801 


207 


10 


489.3900649445 


263 
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3 Other Zeta Functions 



Many other zeta functions which are conjectured to satisfy a Riemann Hypothesis can be ap- 
proximated by sums of hyperbolic gamma functions. Let r(n) denote Ramanujan's r-function. 
Since the function 

oo 

F{z) = T(n) exp(2mnz) (19) 

71=1 

is a modular form of weight 12, we have 

F(ix) = x~ 12 F(i/x), (20) 

which can be used to show 
(2vr)- s r(s) jr ^ = j>2 T ( n ) ( r xS ~ l exp(-2vrnx) dx + f°° x'^x 12 exp(-2irnx) dx\ 

n=l n n=l ^ ^ ' 

(21) 



which implies 



OO / \ oo 

(2^)- 6 -r(6 + iz) J2 53l = E T ^ G ^ 2 ™> 6 )' ( 22 ) 



1 n 

n=l n=l 



again a uniformly convergent sum of hyperbolic gamma functions. (It is known that |t(ti)| = 
0(n 6 ) |Apo90| .) The function defined by the left-hand-side of (|22p is known as the Ramanujan 
H-function, which we denote S/±(z), and the approximate obtained by truncating the series 
on the right-hand-side of ([22]) after N steps we denote Ha,jv(^)- Ramanujan conjectured that 
^a(z) has only real zeros, which is still open. We mention that Ki [Ki0 8] has studied the zeros 
of different approximates to the Ramanujan H-function. 

Conjecture 2 For N G N, ^a,n(z) has monotonic zeros in Q. 

Computations analogous to those described in Remark [T] indicate that the corresponding 
weaker form of Conjecture [2j that there are no nonreal zeros of ^a,n(z) in Q with real part 
smaller than the largest real zero of Sa n(z), is true at least for N < 10. Below we list the 
number of real zeros and the largest real zero of Ha, n (z) for iV < 10. 



K 


largest real zero 


number of real zeros 


1 







2 


9.1937689444922 


1 


3 


13.885647964708 


2 


4 


21.358047646119 


5 


5 


25.047323063922 


6 


6 


28.706422677689 


8 


7 


33.529929734593 


11 


8 


36.535376767485 


12 


9 


40.190608700694 


14 


10 


44.761812314903 


17 
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More generally, we can start with any entire modular cusp form 

oo 

/(z) = £c(n)exp(27riz) (23) 

n=l 

of weight 2k, and set 

c(n) 



B(z) = (2^-^r(k + iz)J2j T L- (24) 

n=l 

The modularity of / can be used to show |Apo90 pp. 137-138] that 

B(z) = J™ f{it) (t k+iz + (-l) k t k - iz ^j dt/t (25) 

OO n QQ J 

YVn) / exp(-u) ({u/2itn) k+iz + {-l) k {u/2Tin) k ~ iz ) — , (26) 

n=l X 7 M 



where we have used the well-known bound c(n) = 0(n k ) to justify the interchange of summation 
and integration. Thus we see that, at least for k even, B{z) is also a uniformly convergent 
(infinite) linear combination of hyperbolic gamma functions. 

We can do a similar calculation for Dirichlet L-series L(s,x) with x a primitive character 
with modulus q. Assume for the moment that x(— 1) = 1, and define 

E(z, X ) = 7r^ 4+iz ^T(l/A + iz/2)g 1 / 4+i2 / 2 L(l/2 + iz, X ). (27) 

From [DavOO, pp. 68-69] we have (using the fact that x( n ) = x(~ n )) 



S(*. X) = E / exp(-n 2 vrt/ (7 ) t 1 /^/ 2 ^ 
n=l ^ * 

°° /"OO 

+ \/9w(x)y]x(n) / exp(-n 2 vrt/g)t 1 
n=l 71 



/A-iz/2 & 



(28) 
(29) 



5C 



^ , r(l/4 + ^/2,vrn7g) 

n=l v ' rL/ 



n=l 
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where yfqw{x) is a certain complex number of modulus one and x( n ) = x( n ) is the conjugate 
character. Furthermore if x(~ 1) = — 1 we get 

^-(3/4+^/2)^3/4 + iz /2) q ^ A + iz / 2 L{l/2 + iz, X ) (31) 



n=l 



/ exp(-n 2 7rt/q)t 3 / A+ ^ 2 -+ (32) 
^(x)E n ^ n )/ exp(-n 2 ^)t 3 / 4 — / 2 ) T (33) 

n=l 

= Vnv(n) r(3/4 + ^ /2 ' 7rra2/g) + (34) 

<v^(x) E »xw r(3 ( ^7 / '! 3 / ?4-w2 /g) ■ 

n=l v ' ' 

where again \iy/qw{x)\ = 1- 

Let F{z) be the function defined by truncating the series on the right-hand-side of (|26|) . or 
both of the series on the right-hand-side of either (|30p or (|34p . after iV steps. If Conjectures 
and [T] and [2] are true, one might suspect that F{z) also has monotonic zeros in Q, although the 
author has not yet done any computations with these more general sums. 

Another interesting question is where the zeros of T(z,a) are, for a a positive real number. 
Neilsen |Nie65j showed that T(z, a) has no zeros in 3t(z) < a, and Gronwall [Grol6| proved that 
T(z,a) has infinitely many zeros in Q. Mahler [Mah30] showed that, as a — > oo, the zeros of 
T(az,a) cluster about the limiting curve 

$l(zlogz + l- z) = 0. (35) 

Tricomi and other authors have investigated the zeros of (the meromorphic continuation of) 
T(z,a) as a function of a, for fixed z. In summary, not much information seems to be known 
about the zeros of T(z, a), for a a fixed positive real number ( although the literature contains a 
number of detailed results on the zeros of the lower incomplete gamma function T(z)—T(z, a)). In 
1998 Gautschi [Gau98] published a nice survey of known results on incomplete gamma functions. 
Computer calculations support the following. 

Conjecture 3 For any fixed positive real number a, the incomplete gamma function T(z, a) has 
monotonic zeros in Q. 

Although some analog of Conjecture [3] may be true for hyperbolic gamma functions, in Section 
[5] we show that there exist some choices of a, b € R, a > for which G(z; a, b) does not have 
monotonic zeros in Q. 

Remark 2 To say a function has monotonic zeros in Q is equivalent to saying that the first 
differences of the imaginary parts of the zeros are all nonnegative. The zeros in Q ofT(z, a), a > 
seem to satisfy the more general property that the k-th differences of the imaginary parts of the 
zeros are positive for k odd and negative for k even, for k < 7 or 8. Thus these zero sets seem 
to have extra structure beyond being monotonic. For linear combinations of hyperbolic gamma 
functions the same phenomena seems to occur for x sufficiently large, which may be due to the 
main term in the asymptotics controlling the zeros. 
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4 Asymptotics 

Throughout this section z = x + iy £ Q, x, y > 0, 8 = arg(z), < 6 < ir/2, a, b £ E, a > 0. We 
begin with Stirling's formula and some other known results: 



^(f)' (1 + 0(1/,)) 



(where — tt < arg(z) < 7r) 



ar 



T(z,a)=T(z)- — (1 + — — + , , . . 9) 
ze a \ z + 1 (z + l)[z + 2) 



+ ...+ 



(« + !)* 



+ 



f |EMUT531 Vol. II, p. 135]) 



T 4±A = Z b (l + 0(l/z)) 



(36) 



(37) 



(38) 



([EMOT53, Vol. I, p. 47]), where in the big-Oh results we mean as \z\ — > oo. From f)36|) we see 
that 




-—exp(xln\z\ -x-y9)(l + 0{l/\z\)) . 



Now T(z, a) = if and only if T(z, a)/T(z) = so by ([37]) T(z, a) = 0, |z| large implies 

a z 



zT(z)e a 



1 



or 



exp(xln \z\ — x — y6)y / 2ir\z\ ~ exp(xlna — a). 



(39) 



(40) 



(41) 



If y remains bounded and x — * oo, the left hand side above grows too fast, while if x remains 
bounded and y — > oo it decreases too fast. Hence we need both x,y — > oo. Furthermore, y0 
must be asymptotic to xlnx, hence 6 — ► it/2 and the zeros of T(z, a) in Q satisfy y ~ -xlnx as 
|z| — ► oo. 

To perform the same analysis for G(z; a, b), first note that 



\T{iz)\ = J-r exp {-y In \z\ + y - z(0 + tt/2)) (1 + 0{l/\z\)) 



\T{-iz) 



z\ 
2tt 




exp (y In \z\ - y + z(0 - tt/2)) (1 + 0(1/|«|)) . 



(42) 
(43) 



Thus using ([38 

|r(6 + »z)| 



~,b+iz I 




z| 6 exp (-y In |z| + y - x(0 + tt/2) - (b - y) hxa) (1 + 0(l/\z\)) (44) 



' U '" ^|z| 6 exp(yln|z| - y + x(0 - vr/2) - (b + y) ma) (1 + 0(l/|z|)) . (45) 
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From (1371) . 



^, ^ TCb + iz) T(b-iz) 

G ^ a > b ) = -^unr 1 + -^ir 1 ( 46 ) 



One finds 



(b + iz)e a £^ Q (b + iz + l) k (6 - iz)e a ^ (6 - iz + l) fc ' 



a K 1 v-^ a 



V V (47) 

(b + iz)e a f^ Q (b + iz + l) k (b-iz)e a ^ (b-iz + l) k 

_2{a-b) 2(6 3 - 2ab 2 - ab + 3a26 + 3a 2 - a 3 ) 6 

If y remains bounded and x doesn't, then the first two terms on the right-hand-side of (|46p 
approach like exp(— xtt/2), so by (|4T|) . the expansion of G(z; a, &) in negative powers of z has a 
nonzero coefficient of z~ 2 , unless a = b in which case it has a nonzero coefficient of z~ 4 . In either 
case it cannot equal for sufficiently large x. If x remains bounded and y doesn't, then the 
second term on the right-hand-side of (|46p blows up, while the others approach 0. So as \z\ —►00, 
for fixed a, b, if G(z; a,b) = we need both x, y — > 00, and thus from (|45p x(# — vr/2) + ?/ln |z| 
cannot approach positive or negative infinity too quickly. In fact the zeros of G(z; a, b) must 
satisfy x ~ ^ylny as |z| — > 00. 

The argument above also applies to any function of the form 

N 



s ^u k G(z;a k ,b k ) (48) 



k=l 

a k ,b k £M,u k £ C,a k > 0, i.e. any C-linear combination of hyperbolic gamma functions. For if 
you have a linear combination of terms like (|46|) . more than one of which is approaching 00, the 
linear combination must also approach 00, since by taking into account the contribution of a k ,b k , 
no two such terms can approach 00 at the same rate. The other parts of the argument follow 
through similarly (the coefficient of l/z 2 i in the appropriate version of (|47p must be nonzero for 
some j, else (j48]) would essentially reduce to a linear combination of Gamma functions, which 
would not be entire) and thus the zeros of any function of the form (|48p also satisfy x ~ ^ylny 
as \z\ — > 00. 



5 Linear Combinations 

The examples in Section [3] involving Dirichlet characters motivate defining, for a,b £ M, a > 0, 
w,a £C, \w\ = 1, the generalized hyperbolic gamma function as 

, \ T(b + iz,a) _T(b-iz,a) , An . 
G(z; a, b, a, w) = a + a ^ V (49) 

The author has made several hundred computer runs, calculating the zeros of various arbitrary 
C-linear combinations of generalized hyperbolic gamma functions in different regions of Q. Sur- 
prisingly, in all of these runs the zeros turned out to be monotonic. Perhaps this results from 
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a mysterious analytic principle, not yet understood, which causes generic sums of generalized 
hyperbolic gamma functions to have a high probability of having monotonic zeros. In this case 
the RH could be a consequence of this analytic principle, combined with Hardy's theorem that 
H(z) has infinitely many real zeros. 

It is not the case though that all linear combinations of generalized hyperbolic gamma func- 
tions have monotonic zeros in Q. For example, consider any function of the form t\A{z)-\-t%B(z), 
where t\,t2 are positive real numbers and A(z), B(z) are any two H-functions, corresponding to 
different zeta functions from Section^ which are real on the real line. Since A(z) and B(z) have 
infinitely many real zeros, as z — > oo along the real line they each oscillate from positive to neg- 
ative. Unless there is an unsuspected correlation between the two, any real linear combination 
of them will also oscillate and thus have infinitely many real zeros. But only for very special 
choices of A,B,t\,t2 will this linear combination correspond to a zeta function with an Euler 
product, and without an underlying Euler product it is generally expected that such functions 
will have infinitely many non-real zeros as well. In particular, 3(2z) + Ea,5(z), the sum of the 
Riemann H- function and the fifth approximate to the Ramanujan 3 function, has a non-real zero 
with real part between 19 and 23, and a real zero at z = 24.99871. (Here S(z) is evaluated at 2z 
to make the two functions compatible, since the expression (|12p of S(z) in terms of hyperbolic 
gamma functions involves z/2, while that of Ha involves z.) 

More simply, one can create an example of non-monotonic zeros by considering what happens 
to Hjv(^) as z — » oo along the positive real line. From (fUJj) and (|4"7|) we get 

G( X -,a,b) = 2 -^^ + 0^y (50) 

Applying this to (|13p yields 

= W( A - - 5/4) J_ + p fl_ N 

exp(n z 7r) x z \ x J 

which shows 

^(x) = -.01974938206/z 2 + 0(l/x 4 ) (52) 
$ 2 (a?) = .01974934121/x 2 + 0(l/x 4 ) 
$ 3 (x) = .4132639753 x lO" 7 /^ + 0(l/x 4 ). 

For k > 3 the coefficient of 1/x 2 in <&k(x) is positive. Thus the coefficient of 1/x 2 in Htv(x) is 
positive for k > 3 and negative for 1 < k < 3. It follows that t&\(x) + &2(x) is positive for 
sufficiently large x if t = and negative for sufficiently large x if t = 1, and so t$i(z) + ^(-z) 
has a zero on the real line for some large x and some < t < 1. One finds in fact that there 
is a real zero between x = 10 5 and x = 10 6 when t = .999997907459, which occurs after many 
non-real zeros, giving an example of non-monotonic zeros. This real zero travels travels quickly 
left as t increases through real values, arriving at x = 39.53248 (the largest real zero of ^(z)) 
when t = 1. 

One could hope that any series of the form 

oo 

F{z) = Ck^k(z), c k e E, c k > (53) 
fc=l 
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has monotonic zeros, modulo the problem of the forced real zero described above, which could 
be avoided by say requiring c\ = C2 = C3, since then 

for some c > 0. A family {fi(z), . . . , fk(z)} of polynomials with real coefficients is called com- 
patible if 

k 

has only real zeros whenever Cj G R, c,- > for 1 < j < k. It is called pairwise compatible if 
{fii fj} forms a compatible family for each pair 1 < i < j < k. Chudnovsky and Seymour [CS07] 
have shown that a family of polynomials whose members have only real zeros and positive leading 
coefficients is compatible if and only if it is pairwise compatible. Does a similar statement hold 
if we replace real polynomials with positive leading coefficients and only real zeros, by even, real, 
entire functions having monotonic zeros in Q? 

Eq. (I50p also leads to an example of a single hyperbolic gamma function with non-monotonic 
zeros. For very large x clearly G(x; a, 2a) < and G(x; a, 0) > 0, so there must be some value 
of 6, < b < 2a for which G(x\ a, b, 1, 1) = for some very large x, which will thus have 
non-monotonic zeros in Q. 

Another interesting phenomena occurs when we consider the zeros of Hfc(z) + t$k+i(z)- As 
we let t vary continuously from to 1, computations indicate that the imaginary parts of the 
non-real zeros decrease monotonically (i.e. continuously), in a very regular manner, until, for 
high enough k, they collide with the corresponding zero (from Schwartz reflection) in the fourth 
quadrant, and arrive on the real line, where they remain. 

Conjecture 4 Fork > 1, the imaginary part of each non-real zero ofz*k{z)+t§k+i{z) decreases 
monotonically (i.e. continuously) as t goes from to 1. 

6 The modulus on vertical rays 

It is known that the RH is equivalent to the statement that the modulus of H(z) is monotone 
increasing along any vertical ray which starts at a point x > on the nonnegative real line and 
travels straight upward to x + ioo. (Clearly if the RH is false the statement is false. On the 
other hand, if the RH is true, start with Hadamard's factorization theorem 

|H(z)| = |S(0)| nil (56) 

i 

where the a« are the positive real zeros of H(z), and take the partial derivative with respect to 
y, where z = x + iy. This is easily seen to be positive for positive y.) 

For a given function F(z), analytic in Q, and a a nonnegative real number, let M(F,a) 
denote the value of y > where the function \F(a + iy)\ is minimal, i.e. the height where the 
minimum of the modulus of F occurs on the vertical ray starting at a and going straight up. 
(If this minimum occurs at more than one y, let M(F, a) denote the lim inf of such y.) Finally 
call the set of pairs {(a, M(F, a))} in Q the "M-curve" of F. 
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Examples of these curves are given in Figures [H [21 and [3] below. To calculate M(F, a) for 
a given pair F and a the author simply calculated the modulus of F(a + iy) at many closely 
spaced grid points y and then chose the y which gave the minimum of these numbers (from 
the asymptotics, the modulus of a sum of hyperbolic gamma functions increases quite rapidly 
when y increases beyond a certain point). This same procedure was followed for several closely 
spaced grid points a, and the pairs (a,M(F,a)) then plotted on a grid using Maple, with the 
result looking something like a continuous curve. They seem to have the (as yet unexplained) 
property that the non-real zeros of F in Q occur at the same places as the local maxima of the 
M-curve. If so, this would give another (rather informal) method for calculating the zeros of 
3jv(^) and other sums of hyperbolic gamma functions which doesn't depend on the argument 
principle or Newton's method. Other computations indicate this property may also hold for any 
real polynomial with no two zeros in the upper half-plane with the same real part. 

6.1 Some Notes on the Computations 

The list of zeros of Ei(z) in the Appendix was calculated using the argument principle and 
Newton's method by a call to the function analytic in Maple, as were the zeros of Ha, 1(2), 
5^,2(2) j and the zeros of the sum of two generalized hyperbolic gamma functions. This procedure 
failed however to compute the zeros of ^(z), as the program never finished even after running 
for over two days on the Sun system at the University of Pennsylvania. Note that Hfc(js) is a 
sum of 4fc incomplete gamma functions, which may explain why the computation of the zeros 
of Ejfe(z) quickly becomes difficult. The author was able to find the zeros of ^(z) by starting 
with the zeros of Hi (2) and using Newton's method to find the zeros of Hi (2) + 1 $2 (2)) for t a 
small positive number, then recursively using these new zeros and Newton's method to find the 
zeros for a slightly larger value of t, slowly increasing t until t equaled 1. The author then tried 
to compute the zeros of ^(z) in the same way, by starting with the zeros of 3z(z) and using 
Newton's method to compute the zeros of B,2(z)-\-t ^s(z) for small t, and gradually increasing t as 
before. This worked well until t became very close to 1, about t = .99, at which point Newton's 
method no longer converged. If the property discussed in the previous paragraph holds though, 
we can see from Figure [2] where the non-real zeros of ^(z) in the range < %l(z) < 120 are. 
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Figure 1: The M-curve (in blue) for ^(z), < $t(z) < 100. The zeros are overlaid in red; note 
they occur at the local maxima of the M-curve. 
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Figure 3: The M-curve for a sum of 8 arbitrarily chosen generalized hyperbolic gamma functions. 
The zeros (of the sum of these 8) are overlaid in red. 
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7 Appendix 



Digits := 25, zeros of first H-approximate Ei(z) in rectangle with corners (0, —.1), (100, 500). 
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04543957882981756479858 










20 


62534600592171760132974 


-1- 


2.697151842339519632505712 


I 


26 


05616693357829946749575 


4- 

i 


7.125359707612690330897455 


I 


31 


50143137824977099308422 


4- 


10 


.72915037105496782822450 


I 


36 


72702276874255239918647 


4- 

1 


13, 


.75961410603683555833019 


I 


41 


73703479849622101486046 


4- 


16 


.44012737324329251859479 


I 


46 


56622866997881255099908 


4- 


18, 


.88186965378958902053812 


I 


51 


94458589^1 1 fi9Q45S4fi8QQ0 


j_ 


21, 


. 14750420601374895347492 


T 

_L 


55 


.79525368022472028456165 


+ 


23 


,27625685820891335493023 


I 


60 


.23621426525993802296865 


+ 


25 


.29458549895993860216014 


I 


64 


,58150497097301796850798 


+ 


27, 


.22133555778112035831075 


I 


68 


,84235653395121330563843 


+ 


29 


.07049609150585601287785 


I 


73 


.02789933182939276748060 


+ 


30 


.85279227139366634464017 


I 


77 


, 14567324003250763696303 


+ 


32 


,57666324204392832752644 


I 


81 


,20199121212953110713480 


+ 


34. 


.24889253114939152723783 


I 


85 


.20220345212231662722890 


+ 


35. 


.87503096670553315342957 


I 


89 


,15089297349449318064800 


+ 


37 


.45968995259880236581690 


I 


93 


,05202292717284600209187 


+ 


39. 


.00675061925213970478000 


I 


96 


,90904939663401491219210 


+ 


40 


.51951660155401879741380 


I 



1st differences of imaginary parts of zeros 

2.697151842339519632505712 

4.428207865273170698391743 

3.603790663442277497327045 

3.03046373498186773010569 

2.68051326720645696026460 

2.44174228054629650194333 

2.26563455222415993293680 

2.12875265219516440145531 

2.01832864075102524722991 

1.92675005882118175615061 

1.84916053372473565456710 

1.78229617988781033176232 

1.72387097065026198288627 

1.67222928910546319971139 

1.62613843555614162619174 

1.58465898589326921238733 

1.54706066665333733896310 

1.51276598230187909263380 

Digits : = 35, zeros of £2(2) i n rectangle with corners (0, —.1), (102, 100). 

14.1347251016150223590867934323428320 
21.0220425550989420016995644399118565 
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25.0108186655252580507941448128443705 
30.4267684045343046403275534114141879 
32.9244008910878391337889771099093884 
37.8603410698702633339745855804483258 
39.5324810797879508571153350051479630 

43.1389080680988950956236929709951507 + 3.28097100306881350163884998294539870 I 
47.5227563037011494043923449534796642 + 6.25250911174970354129046143001000600 I 
51.8283147335875594556295444210165891 + 8.95857384762622450769949019637180856 I 
56.1120029199774070873341435705058778 + 11.4796085175258061625014962442075566 I 
60.3502919249623580499827793888173447 + 13.8416168548275384327656057308155260 I 
64.5390675686875954205140494065924195 + 16.0705629583418306170962131727661418 I 
68.6761126675575950561373116041091163 + 18.1868197094425056650011185673044960 I 
72.7617836933667997393215687715081074 + 20.2066697023948159796720721861483021 I 
76.7975312791721943222372025862518267 + 22.1432248556299027619859319905393103 I 
80.7854051988095067379835000872591797 + 24.0071081808438939881307280663316996 I 
84.7277019646143359935370107950650673 + 25.8070157111458796320542608491379890 I 
88.6267772741056992788936274593068783 + 27.5501321650503597200806915781312903 I 
92.4849402660525200867715711075727066 + 29.2424518834035459512085274002314648 I 
96.3043979005927443711172685790765488 + 30.8890225459009695084228517057586202 I 
100.087228624602886284768301012877784 + 32.4941324878154139135000560357642324 I 

Digits := 35, zeros of £3(2) in rectangle with corners (0, —.1), (66,67). 

14.13472514173469376946955013374227983392383419451107438 
21.02203963877155659023693964971327708279243073624430253 
25.01085758014566673273814860967472222032345345332059477 
30.42487612586063748020316156144982639562443306193294638 
32.93506158773265354499645423632051649764195395468931578 
37.58617815896445021695716617231461056957864550061946742 
40.91871901006637538924149719014594063254432154184503222 
43.32707329088964789958258005109768490700030014005496697 
48.00515051035397786165873118175586051353275737166717210 
49.77383399904302558365984997979234680103752864503720614 
52.97031190983202078103080161894128897412215226172440818 
56.44637829985894482385431963634767652784719134486888209 
59.34511184966886217051112646891854833040742035322093409 
60.83672336038805375320166674337805885068490701004798353 
65.03207377198913910137679883482704065081967504834437506 

Digits := 35, zeros of first approximate Sa,i to the Ramanujan H-function in rectangle with 
corners (0, -.1), (100, 100). 

9.5489635091412543125672429455388951 + 1.7119172216212706912355445109132207 I 
13.789433051893611052566335275629933 + 5.2848734079809646491896114233647142 I 
17.427358570320731571246782104695096 + 7.7540563400814520087583759352815352 I 
20.757098189027862070211629799284288 + 9.7974743824676388131629735365802915 I 
23.879056232416349736892490077064465 + 11.592456287772395127573044664358961 I 
26.847356442625470267126707245453301 + 13.219199923922411333605849186386965 I 
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29.695457505212671889597523331405240 + 14.722240962151850502838429485181816 I 
32.445858241734786458856784110868772 + 16.129361168830359821304709531815352 I 
35.114596915766277030958889001310311 + 17.459270918804088444729979963005156 I 
37.713610437352047890893961093475018 + 18.725280350575037459960230440611499 I 
40.252083817737620602530239123250025 + 19.937256671472885578681886805140642 I 
42.737274291251073747956340698712113 + 21.102754510914081337479910411119165 I 
45.175041128226087752248534272004155 + 22.227709943351466358153291572887016 I 
47.570200654826011921382328419781717 + 23.316887844914435394477534028073124 I 
49.926772425743961667810425669964510 + 24.374181909750589494277749519061023 I 
52.248154899971609917404551354946461 + 25.402822596146688919661995414148711 I 
54.537253919041181821883861911330465 + 26.405525330914909309258552834547793 I 
56.796578677359823555284912818063663 + 27.384598689419083619944995509463388 I 
59.028314745351122620881979673545968 + 28.342025013215928897705357026998109 I 
61.234380542300284878652304849348298 + 29.279521587539252797178004730959225 I 
63.416471643946175987038219068629433 + 30.198587816232713426722974287579177 I 
65.576095995892802432660080685311982 + 31.100542121402673258701892990637098 I 
67.714602225166643613493231218811692 + 31.986551176775085187674532167087986 I 
69.833202642000750529879098123320031 + 32.857653335480587110184136136077900 I 
71.932992106046633644353345531169636 + 33.714777601864923783681632805009935 I 
74.014963635252309537091144416299462 + 34.558759141214355169362183745447748 I 
76.080021422732136425762615907263059 + 35.390352069528146990609600601997921 I 
78.128991771589170231983397657306861 + 36.210240084528180967081058119582790 I 
80.162632342789627094434403795275856 + 37.019045367211828466413017905373453 I 
82.181640025246587241041487951836980 + 37.817336085881271861478883874617446 I 
84.186657672248860977039557949370316 + 38.605632761831562562091441356062965 I 
86.178279898670686132064566105956960 + 39.384413700925457309936759214348925 I 
88.157058095044852280187866298051226 + 40.154119653347769812443058137516104 I 
90.123504784721498632008498428129335 + 40.915157831527436633928427478875093 I 
92.078097426893282987651528394221783 + 41.667905391109706656241503132525470 I 
94.021281749721926032933713327889188 + 42.412712460188349826344885668577328 I 
95.953474683021781433616191339406534 + 43.149904786473380482490406058490112 I 
97.875066948097497309259057547588008 + 43.879786059713960896674731680251557 I 
99.786425352756449975929052839495573 + 44.602639956801403127156000977181747 I 

Digits := 25, zeros of second approximate Ha, at to the Ramanujan H-function in rectangle 
with corners (0, - .1), (100, 200). 

9.19376894449217893 

13.8348516743546257 + 1.54300060155520050 I 
17.0006907368319637 + 4.58047657550007841 I 
20.1893418670461022 + 6.99896197324829350 I 
23.2665587059084961 + 9.04437482423736081 I 
26.2157755564237380 + 10.8464634984827865 I 
29.0559366154717340 + 12.4806210190923921 I 
31.8034718883051467 + 13.9910040730670684 I 
34.4718383897903041 + 15.4050720323654861 I 
37.0718504676122273 + 16.7414677662025767 I 
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39 


6120959833238977 


+ 


18 


0135340186718953 


I 


42 


0994917987503797 


+ 


19 


2311725419877171 


I 


44 


5397026159453686 


+ 


20 


4019721238147549 


I 


46 


9374285351558464 


+ 


21 


5319037938006434 


I 


49 


2966167316589794 


+ 


22 


6257645644359874 


I 


51 


6206201858570109 


+ 


23 


6874764176583670 


I 


53 


9123161924908609 


+ 


24 


7202947836257462 


I 


56 


1741960636482404 


+ 


25 


7269571656683657 


I 


58 


4084343572848694 


+ 


26 


7097915965161903 


I 


60 


6169430187982438 


+ 


27 


6707975724897949 


I 


62 


8014141732518389 


+ 


28 


6117075511920988 


I 


64 


9633543292524033 


+ 


29 


5340343986557084 


I 


67 


1041120187127688 


+ 


30 


4391085162350689 


I 


69 


2249003477575950 


+ 


31 


3281072698702158 


I 


71 


3268155548021774 


+ 


32 


2020785869307989 


I 


73 


4108524066057312 


+ 


33 


0619600720511948 


I 


75 


4779170679704262 


+ 


33 


9085946393025650 


I 


77 


5288379348185691 


+ 


34 


7427434066455078 


I 


79 


5643748117146020 


+ 


35 


5650964167873912 


I 


81 


5852267334967588 


+ 


36 


3762816159807758 


I 


83 


5920386687224333 


+ 


37 


1768724246086591 


I 


85 


5854072948756678 


+ 


37 


9673941604083026 


I 


87 


5658859982298918 


+ 


38 


7483295199556847 


I 


89 


5339892223226212 


+ 


39 


5201232818472100 


I 


91 


4901962662154965 


+ 


40 


2831863625189748 


I 


93 


4349546156253617 


+ 


41 


0378993303884662 


I 


95 


3686828755524769 


+ 


41 


7846154642047194 


I 


97 


2917733614022269 


+ 


42 


5236634258506564 


I 


99 


2045943961865493 


+ 


43 


2553496053959862 


I 



Digits := 25, zeros, in rectangle with corners (10000.0, 500), (10025.8, 10000), of linear 
bination of two generalized hyperbolic gamma functions 

2 

,B k ,otk,w) 

k=i 

where 

0i = 1.0 + .3i, = —3 — i, a.\ = 1.0 + .3i, a 2 = —3 — i, 
w = 0.7648421872844884262558600 + 0.6442176872376910536726144/ 
B x = 1.0, B 2 = 2.2, A 1 =tt,A 2 = 7.853981633974483096156608 

10000.08899697917243762827 + 1943.699378898319010297001 I, 
10000.84424695361683674126 + 1943.828178641310676925944 I, 
10001.59949024766823939791 + 1943.956976210289663002380 I, 
10002.35472686200573704240 + 1944.085771605493580928470 I, 
10003.10995679730829337130 + 1944.214564827159996509445 I, 
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10003.86518005425474436858 + 1944.343355875526428966768 I 
10004.62039663352379834061 + 1944.472144750830350951281 I 
10005.37560653579403595119 + 1944.600931453309188556355 I 
10006.13080976174391025650 + 1944.729715983200321331031 I 
10006.88600631205174674006 + 1944.858498340741082293155 I 
10007.64119618739574334772 + 1944.987278526168757942515 I 
10008.39637938845397052261 + 1945.116056539720588273964 I 
10009.15155591590437124002 + 1945.244832381633766790543 I 
10009.90672577042476104237 + 1945.373606052145440516603 I 
10010.66188895269282807412 + 1945.502377551492710010906 I 
10011.41704546338613311658 + 1945.631146879912629379767 I 
10012.17219530318210962295 + 1945.759914037642206290096 I 
10012.92733847275806375302 + 1945.888679024918401982558 I 
10013.68247497279117440814 + 1946.017441841978131284633 I 
10014.43760480395849326601 + 1946.146202489058262623716 I 
10015.19272796693694481556 + 1946.274960966395618040195 I 
10015.94784446240332639171 + 1946.403717274226973200536 I 
10016.70295429103430821021 + 1946.532471412789057410354 I 
10017.45805745350643340246 + 1946.661223382318553627526 I 
10018.21315395049611805021 + 1946.789973183052098475038 I 
10018.96824378267965122043 + 1946.918720815226282254485 I 
10019.72332695073319499999 + 1947.047466279077648958683 I 
10020.47840345533278453047 + 1947.176209574842696284945 I 
10021.23347329715432804280 + 1947.304950702757875648075 I 
10021.98853647687360689210 + 1947.433689663059592193412 I 
10022.74359299516627559224 + 1947.562426455984204809882 I 
10023.49864285270786185070 + 1947.691161081768026142988 I 
10024.25368605017376660312 + 1947.819893540647322607883 I 
10025.00872258823926404804 + 1947.948623832858314402362 I 
10025.76375246757950168151 + 1948.077351958637175519890 I 

1st differences of imaginary parts 

0. 128799742991666628943 
0.128797568978986076436 
0. 128795395203917926090 
0.128793221666415580975 
0. 128791048366432457323 
0.128788875303921984513 
0.128786702478837605074 
0.128784529891132774676 
0.128782357540760962124 
0.128780185427675649360 
0.128778013551830331449 
0.128775841913178516579 
0.128773670511673726060 
0. 128771499347269494303 
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0.128769328419919368861 
0.128767157729576910329 
0. 128764987276195692462 
0.128762817059729302075 
0. 128760647080131339083 
0.128758477337355416479 
0.128756307831355160341 
0. 128754138562084209818 
0. 128751969529496217172 
0.128749800733544847512 
0.128747632174183779447 
0.128745463851366704198 
0. 128743295765047326262 
0.128741127915179363130 
0.128738960301716545337 
0.128736792924612616470 
0. 128734625783821333106 
0. 128732458879296464895 
0.128730292210991794479 
0.128728125778861117528 

2nd differences of imaginary parts of zeros (times minus 

0.2174012680552507 x 10~ 5 
0.2173775068150346 x 10~ 5 
0.2173537502345115 x 10~ 5 
0.2173299983123652 x 10~ 5 
0.2173062510472810 x 10" 5 
0.2172825084379439 x 10~ 5 
0.2172587704830398 x 10~ 5 
0.2172350371812552 x 10~ 5 
0.2172113085312764 x 10~ 5 
0.2171875845317911 x 10" 5 
0.2171638651814870 x 10~ 5 
0.2171401504790519 x 10~ 5 
0.2171164404231757 x 10~ 5 
0.2170927350125442 x 10~ 5 
0.2170690342458532 x 10~ 5 
0.2170453381217867 x 10~ 5 
0.2170216466390387 x 10~ 5 
0.2169979597962992 x 10~ 5 
0.2169742775922604 x 10~ 5 
0.2169506000256138 x 10~ 5 
0.2169269270950523 x 10~ 5 
0.2169032587992646 x 10~ 5 
0.2168795951369660 x 10~ 5 
0.2168559361068065 x 10~ 5 
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0.2168322817075249 x 10~ 5 
0.2168086319377936 x 10~ 5 
0.2167849867963132 x 10~ 5 
0.2167613462817793 x 10~ 5 
0.2167377103928867 x 10~ 5 
0.2167140791283364 x 10~ 5 
0.2166904524868211 x 10~ 5 
0.2166668304670416 x 10~ 5 
0.2166432130676951 x 10~ 5 
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